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ABSTRACT 

Faraday rotation (FR) is widely used to infer the orientation and strength of mag- 
netic fields in astrophysical plasmas. Although the absence of electron-positron pairs 
is a plausible assumption in many astrophysical environments, the magnetospheres of 
pulsars and black holes and their associated jets may involve a significant pair plasma 
fraction. This motivates being mindful of the effect of positrons on FR. Here we derive 
and interpret exact expressions of FR for a neutral plasma of arbitrary composition. 
We focus on electron-ion-positron plasmas in which charge neutrality is maintained 
by an arbitrary combination of ions and positrons. Because a pure electron-positron 
plasma has zero FR, the greater the fraction of positrons the higher the field strength 
required to account for the same FR. We first obtain general formulae and then specif- 
ically consider parameters relevant to active galctic nuclei (AGN) jets to illustrate the 
significant differences in field strengths that FR measurements from radio frequency 
measurements. Complementarily, using galaxy cluster core plasmas as examples, we 
discuss how plasma composition can be constrained if independent measurements of 
the field strength and number density are available and combined with FR. 

Key words: magnetic fields; plasmas; galaxies: jets; radio continuum: general; ISM: 
magnetic fields; methods: observational 



1 INTRODUCTION 

Faraday rotation (FR) describes the effect by which the elec- 
tric field vector of propagating electromagnetic (EM) radia- 
tion rotates as it propagates through a magnetized plasma. 
The EM waves interact with charged particles such that the 
left and right handed components experience different re- 
fractive indices and propagate at different phase velocities. 
The net electric field vector rotates as the wave propagates. 
The total amount of this FR between source and observer 
depends on the strength and orientation of the interven- 
ing magnetic field along the line of sight, the plasma den- 
sity and the plasma composition. FR has been used, for 
example, to constrain the the stre ngth, gradient, and direc- 
tion of m agnetic field i n AGN iet s('Asada et al]l200l [iooi : 
lAttridge e t al. 2005; Zavala fc Tavlor 2005). galaxy clusters 
(REF) and the Galactic ISM (REF). 

Typically, calculations of FR are made for a neutral 
ion-electron plasma. While this is a good approximation for 
the plasma of the Galactic ISM and galaxy clusters, it is 
less reliable for magnetospheres and outfiows around neu- 
tron stars and black holes. In particular, constraining the 
composition of jets in active galactic nuclei (AGN) is obser- 
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vationally difficult because only the electrons or positrons 
radiate efficiently so free energy in protons is hard to di- 
rectly to detect. The plausible predominance of electron- 
plasma in the magnetosphere of the black hole and the un- 
certainty as to whether jets emanate directly from its magne- 
tosphere (|Blandford fc Znaiekl 1977|) or from the surround- 
ing accretion disk jLivio. Qgilvie fc Pringlelll999h has con- 
tributed to making the composition determination a long 
standing puzz le. The b asic question of how far a jet would 
be force-free (|Li et al ] 12006)) is also a fundamental issue. 
Indirect theoretical constraints o n AGN jet plasm a com - 
positi on reach mixed conclusions: ICelotti fc Fabiaiil (llQQST ) 
and Ghisollini (200 ^ favor predominately electron-proton 
jets, whilst Revnolds et al.l (fl996i) favors a predom i nant! 
electron-positron jet for M87. Hubbard fc BlackmanI (|200i 
argued that on the largest scales, stellar wind mass loading 
will significantly proton load the jet regardless of its initial 
composition. It is therefore important to identify possible 
new techniques for determining jet plasma composition and 
to be aware of its effect on measurements such as FR. 

With the above motivations in mind, we calculate the 
effect of plasma composition on FR from first principles. In 
section 2 we first derive the general exact formula for FR in 
an arbitrary neutral plasmas. We then consider the special 
cases of a pure ion electron plasmas and an ion-electron- 
positron plasma. In section 3 we solve the exact and high 
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frequency approximations to the general FR equation and 
show quantitatively how plasma composition leads to de- 
generacies in the electron number density rie, the magnetic 
field B. We discuss how these degeneracies can in principle 
be used to constrain the plasma composition in AGN jets, 
and conclude in section 4. 



2 GENERALIZED FARADAY ROTATION 

2.1 Formalism for Arbitrary Neutral Plasmas 

To formally derive FR for an arbitrary neutral plasma, we 
assume a cold neutral plasma in a background external mag- 
netic field Bex , subject to a perturbation from a propagating 
electromagnetic wave. (The cold plasma approximation for 
FR has been shown to be effective for the ele ctron contribu - 
tion to FR even for quasi-relativistic plasmas (jSkillindllQTlT l. 
and we discuss this further below Ea. (|16p ). For the electric 
field E, magnetic field B and induced particle velocity Vs 
(where the index s indicates particle species) we write 

E = Ei, 

B = Bl + Bex, 

Vs = Vsl, (1) 

where Ei and Bi are perturbations such that \Bi/Bex\ ^ 
1, \E-i/Bex\ <C 1 and Vsi\ <C c. We also assume a neutral 
plasma so that 



UsOSs = 0, 



(2) 



where nso is the unperturbed density of a particle of species 
s, and Bs is charge of particle of species s. 

Using the above formalism . Max well equations become 
(e.g. iGurnett fc Bhattacharie3 (|2005l ll: 



V- (3) 



„ - IdBi „ in ^ IdEi 

yxEi = --, VxBi = — J-f- — 

c ot c c at 

If the current density J = X/ ^aoe-s^si = 0, then Ei & 
B\ are decoupled, resulting in the plane wave vacuum equa- 
tions. However a finite J and the Lorentz force equation 



dvs , ^ 1 ^ ^ 

ms—— = CsiEl + -Vs X Be 

at c 



(4) 



imply that in general, that E, B, and Vs are all coupled. 

To quantify the interaction between the particles and 
EM fields, we take Vs,Ei oc e'^*'^"'^*) so Eq.Q becomes 



■^{E, + ^-XBex) 
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(5) 



where uj^s = esBex/msC is the cyclotron frequency of species 
s. The current density can then be expressed as the product 
of a conductivity tensor and the Ei field, n amely J = if -Ei, 
where the conductivity tensor is given by (lAshokll2004 ) 



msLo\i-(^Y] y 



~ ^ijk^cs,k 



—iurUsVsy = es{Eiy —B), 

c 

—iujmsVs^ — esEiz, 
and, the conductivity tensor becomes 



(6) 



We now take Bex = (0, 0, B) so that the components of 
Eq.([5ll become 



-iumsVsx = es{Eix + —^B), 
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Eq.(13| with \7 ^ ik then becomes 
■■2 47r^, 
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El = 0. 



(7) 



(8) 



(9) 



from which the secular equation for the FR effect is 

-iD \ / Eix \ 
S~n^ Ely =0, (10) 

P J \ El, J 

where 5 = 1 - Z^^' ^ ^^^^^y P^^- 
and the plasma frequency and wave vector are given 
respectively by LUps — ^^"""^^ and k = ^ z. 




The solution of Eq. (|T0|) for refractive index n produces 
non-trivial FR when the the two solutions for n are dis- 
tinct, corresponding to left and right handed polarizations; 
n\=S — D, njf=S + D with two associated Ei field eigen- 
vectors. The transverse {x,y) components of E for each re- 
fractive index are of the same magnitude but have differ- 
ent phases, that is. El = (Eq, —iEo,0), Er = (_Eo, iiSo, 0), 
where ±i arises from the differentiation of velocity and posi- 
tion over time in the Lorentz force law. The equal amplitude 
of the transverse E field components then imply circularly 
polarized waves. 

The different phase velocities {c/tll, c/nB.) cause the 
propagating left and right handed circularly polarized waves 
to experience a net phase angle difference when they prop- 
agate over the same distance. As a result, the net electric 
field phase angle (0 — Tan~^ {Ey / Ex)) that comes from the 
superposition of these handed waves rotates along the prop- 
agating distance. This is the FR. The change (j) along the 
propagation distance is 



dz 




1 - 



E 



Having derived the general formalism for a neutral 
plasma of arbitrary composit ion and the exact equ ation for 
FR (Eq.dlT}), we note that iHall fc Shuklal (|2005l l consid- 
ered FR in an ion-electron positron plasma producing the 
approximate analytical result 



d4 
dz 



2ne-^B 



(12) 



where Zi is the ion charge number. Ea. p2|l agrees with 
Eq. (|ll|l in the high frequency limit for an ion-electron- 
positron plasma, a point we will return to in section 2.3. 
Ea. (|12p indicates there is no FR in case of an electron- 
positron pair plasma(ni = 0). For Ui 7^ 0, ions generate 
the FR by breaking the symmetry of a pair plasma. 



Effect of Plasma Composition on the Interpretation of Faraday Rotation 3 



2.2 Ion-electron plasma 

For a pure (hydrogen) ion-electron plasma, Eq.([2ll takes the 
form Ue + Ui = 0, where Ue and rii are the electron and 
hydrogen ion number densities. The summation over s in 
Eq. (|lH) also involves terms corresponding to electrons and 
ions respectively. However, because of the large ion to elec- 
tron mass ratio rwi/me = 1836, the io n terms are typi- 
cally i gnored (e.g lAsada et al.l l|2002l . [20o3 ): IZavala fc Tavlorl 
Then the rotated angle integrated along the line of 
sight for the distance / (e.g. Rcvnolds et al. (1996)) becomes 



2c VV Lj{lU — LJce) 



1 - 



Tie B cos 9 dz 



V 27rmic4 7o 



n^B cose dz] ■ y- = RM ■ 



■ dz 



(13) 



where the second relation follows for lj S> oJpc , t^ce ■ The 
general procedure for determining RM is to measure at 
multiple wavelengths and infer a slope of the cf) vs. A line. 



2.3 lon-electron-positron plasma 

For the hydrogen ion-electron-positron case Eq.© becomes 

— eue + en^+ + em — 0, (14) 

where nj.+ is the positron number density. We now define 
X = /we so that 

n^+^UeX, Hi ^ ne{l - /Ue) ^ ne{l — X). (15) 
Eq. (|lip then becomes 



dz 2c 




ljj{u) — U)ce) Ll!{uJ — ce+ ) l^iy - 



plasma term. However, for a pair plasma in which the elec- 
tron and positron cold plasma terms cancel exactly, warm 
plasma correction terms would not cancel exactly and a fi- 
nite contribution would remain as the positron and electron 
correction terms do not cancel. We ignore these small cor- 
rections for present purposes and leave further discussion for 
future work. 

In the high frequency limit, qi, ^ 1 and qn ^ 1, and 
we can approximate Ea. (|16|l as 



dz 



■n,B{l- X) 



= B(1-X)(-^ 



( — 

^ TYIpC ' 



183621.^2 



18362l 



(17) 




where qL represents the second term under the first square 
root of the second equality and g_R represents the second 
term under the second square root of the second equality. 

Ea. (|ll|l . (I12|) and (|16[) presume a cold plasma and 
it is instructive to comment on the validity of these ex- 
pressions for a warm plasma. For the latter, motions of 
charged particles are influenced by thermal effects in addi- 
tion to the electromagnetic force and to express the current 
density the solution of Vlasov equation is necessary (e.g . 
iGurnett fc Bhattacharied (|2005l ). chapt. 9). ISkillingj (Il97lh 
studied FR for a warm ion-electron plasma and found that 
for large frequencies away from resonances, the correction to 
the electron contribution to FR is small compared to the cold 



Using Eq.dlSj, then Eq.([T7| is the same as Ea. (fT2)) . 

For a pure neutral pair plasma (X=l), the right side 
of Ea. (|16p (or (1171) ) vanishes. The FR vanishes because the 
equal mass of positrons and electrons induce the same phase 
speeds for oppositely handed EM waves. This contrasts the 
limit of the previous subsection of a pure neutral (hydrogen) 
ion-electron plasma {X = 0), for which the mass asymme- 
try leads to unequal phase speeds of the oppositely handed 
waves and a flnite right side of Ea. (|16|l . In general, for 
^ X 1 with rie , -Bea; and source distance fixed, the 
right side of Ea. (ll6p decreases with increasing X. We dis- 
cuss solutions of Ea. p6|l in the next section. 

Note that the exact FR expression (I16|l has singularities 
when the wave frequency of the EM wave coincides with the 
particle cyclotron frequencies, i.e. at = eB /m^c = 1.76 x 
B (for electrons and positrons) and eB/niiC = 9571 B 
(for ions). The FR would exhibit sharp resonance features 
near the singular points, allowing the B field to be inferred 
in principle. However, for applications to extended jets of 
AGN and larger scale systems, these resonant frequencies are 
generally small compared to the relevant ~ GHz frequencies. 



3 ASTROPHYSICAL IMPLICATIONS 

3.1 General Implications of Plasma Composition 

Vigs.l and 2 show solutions to the exact expression Eq. p6|l 
And the approximation Eq.((T7} for rie, B and X at fixed 
values of RMs (2592 rad/m^ in Fig.l. and 1000 rad/m^ in 
Fig. 2). The RMs were converted into rotation angles for cm 
scale at A= 1.35 cm(22.2GHz). The FR in Fig.l corresponds 
to 3.63 pc (749 Mpc x 0. 5 mas , CI region) from the core of 
6fGN jet 3C 273 f rom IZavala fc Tay l or (2 0^. From syn- 
chrotron emission, ISavolainen et al.l (j2006l ) calculated the 
total magnetic field to be {B ^ 0.06 G) in this region. If 
this were the line of sight field. Fig 1(b) shows that rie ^0.05 
cm~^ for RM (2592 rad m"^)^ the minimum occurring at 
X = Q. Fig 2(c) shows a complementary example for values 



appropriate for a typical X-ray cluster (jjaffel (|l980l )) with 
rie ~ 0.003 cm"^ and RM ~ 1000 rad m"^. The standard 
assumption that X = for a kno wn distance leads directly 
to the inference that B ~ l»G (|Jaffel Il980l ). But for any 



X < 1, Fig 2(b) for example, shows how much stronger the 
field could be. 



For fixed values of RM, Figs 1(a) 



and 2(a) show that 
We and B behave oppositely for each value of X: As in- 
creases (decreases), B decreases (increases). These trends 



4 Kiwan Park and E. G. Blackman 






(a) log B VS Tie 



(b) log B VS X 



(c) log Tie VS X 



Figure 1. Plots of rie, B and X. All panels use a fixed RM = 2592 rad m-2 in Zavala & Taylor (2001) {Fig.3). The line of sight B 
field is estimated at ~ 0.5 mas (~ 3.63 pc) from the core. All of the plots were made using both Ea. (|16| l and Eq.^TT} with A=1.35 cm 
(22.2 GHz). The plots of each equation were overlapped highlighting the efficacy of the approximate equation. The original observation 
frequency is 8GHz but any frequency gives the same RM. 





(a) log B VS Jie 



(b) log B VS X 



(c) log Tie VS X 



Figure 2. Same as Fig.l. but using RM~1000 rad m ^, path length I =300 kpc, in an X-ray cluster core Jaffe (1980) with A=1.35 cm 
(22.2 GHz) Again both Eq.((T6]l and Eq.lfTTll were used for the plots. 



reflect that the RM is roughly proportional to the field and 



the density. Figs 1(a) 1(b) 2(a) and 2(b) show that as X 
increases (decreases), B and increase (decrease) respec- 
tively. These trends result because an increasing X means 
a higher fraction of pair plasma. The latter contributes zero 
FR so that a higher B or is needed for a fixed RM. 



Figs 1(c) and 2(c) reflect these same trends. 

In all panels, the lines resulting from Eq. (|16|l are in- 
distinguishable from those obtained using Ea. (ll7|l for the 
parameters used, highlighting the efficacy of the latter. 



Figs 2(a) 2(b) and 2(c) are very similar to Figs 1(a) 1(b) 



and 1(c) as only the vertical axis scales are different due to 



the different RM choices. Overall, the figures show how an 
unknown plasma composition X implies degeneracies in the 
values Ue and B, or complementarily, how independent mea- 
sures of rie and B can be combined with an RM to constrain 
X. 



3.2 Further Discussion of Applications 

If the RM and two of the three quantities rie, X and B are in- 
dependently known, the FR equation is exactly determined. 
However, even if only one variable and RM are known, the 
other variables can be constrained. 

There have been efforts to interpret the RM with a sub- 
set of in dependently measured v ariables rie, B and T. For 
example, IZavala fc Tavloij l|2002l ) calculated the B field in 
M87 using the RM and an independently determined Ue- 



For RM = -4000 rad m~^, Eq.((T3I was used to get S ~15 
^G, but X = Q was assumed. We can revisit the inter- 
pretation of the measured RM without assuming X = Q 

the cross point of B 



15 ^lG 
^ (dashed 



a priori. As seen in Fig 3(a) 

and RM = -4000 rad m"^ with rie = llOOcm" 
line) is X = 0. In contrast, if we suppose that B = 200 
^lG, which corresponds to the thermal equipartition condi- 
tion when T ~ 10* if {uekBT = B^/Stt where B is the 
moan line of sight field that causes FR, e.g. iGabuzda et al.l 
(2008)), different X values resuh: For RM = 9000 rad m"^ 
and -4000 rad m~^, the X values are 0.83 and 0.92 re- 
spectively. The corresponding values are 913 and 1012 
cm~^. If instead magnetic pressure dominates (e.g. domi- 



nant {j3 = Pg/Ps = 87rn£ 
shown in Fig.3. 



ksT/B <1), then X increases as 



As another set of examples showing the degeneracy 
between B and X, we consider the independently mea- 
sured RM, He, B arid T for clusters A400, A 1795 and 
A2199 lEilek fc Owenll2002l ) which respectively include ra 



dio sources 3C75, 4C46.42 and 3C338 to obtain Fig 3(b) 



B is also the mean line of sight field but we note that the 
definition of the magneti c pressure (Pb = 3 < B \\ >^ /Stt) 
is different from that of I Zavala fc Tavlo j l|2002l ). We have 
taken this into account when interpreting their respective 
data. There, instead of assuming X = 0, as is normally done 
to obtain the B, we chose selected field strengths (straight 
lines on the plot) and identify the constraints this places on 
X by where these lines intersect with the curves. For ex- 
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X X 
(a) log B vs X for M87 (Zavala & Taylor 2002) (b) log B vs X for Abell clusters (Eilek & Owen 2002) 

Figure 3. These LogB vs X plots were made using RM, B, and T from Zavala & Taylor (2002) and Eilek & Owen (2002). The 
rotation angle for (a) was converted from the RM -4000 rad m^^ (dashed line) and 9000 rad m^^ (solid line) at A= 1.35 cm (22.2 GIfz). 
For (b), each data point for RM (50, 1500, 750 rad in'^); rie (0.0021, 0.064, 0.02 cm-^); and T (1.5, 5.1, 1.2 keV) where the values in 
parentheses are for A400, A1795, A2199 respectively, were also plotted for 22.2GHz. The values of /3 correspond to the particular cluster 
which the spearate lines of B intersect, explaining why different values of the field correspond to the same j3 = 0.1 for A400 and A2199. 



ample, for A400, if magnetic and thermal pressures were in 
equipartition (/3 = 1), the data would imply X = 0.55 (Here 
Tie = 0.0021 cm"^ = 1.155 x 10"^ cm"^ m = 9.45x 
10~* cm~^)). However, if magnetic pressure were domi- 
nant (e.g. /3 = 0.1), then X = 0.86 (rie = 0.0021 cm"^, 
= 1.806 X 10"^ cm"^ m = 2.94x 10"" cm~^). The 
extent to which X differs form zero in clusters could depend 
on how close to a radio source core the RM is measured as 
one would indeed expect X = far away from radio jets. 
As better spatial coverage of both FR and synchrotron con- 
straints become available for jets, the use of plots such as 
those of Fig.l and 2 become more powerful. 



4 CONCLUSIONS 

Assuming plasma neutrality, we have generalized the cal- 
culation of FR for a neutral plasma containing electrons, 
positrons, and ions. Positrons contribute equally and op- 
positely to electrons under the approximations considered. 
Thus as the ratio of positrons to electrons X increases from 
to 1, the FR weakens for a given line of sight magnetic field 
and density. Correspondingly, for a given RM, a larger value 
of X would imply a field strength and or electron density in 
excess of that inferred for X = 0. 

While = is often assumed for astrophysical plas- 
mas, this assumption is not necessarily valid for jets and 
magnetospheres of black holes and neutron stars, where in 
fact it is of interest to independently determine X. We have 
shown quantitatively the degeneracies in the space of pa- 
rameters rie, X, and B and how they can be constrained. In 
principle, if independent data on B and ne can be obtained, 
then a given RM measurement can be used to obtain X. The 
pursuit of X has been particularly elusive in the context of 
AGN jets and we hope that the calculations herein provide 
a tool toward this goal, and help motivate the pursuit of 
further data. 
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